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Abstract 

We show that long distance contributions to the rare decays K ttvv and K — ^ vrf+f" 
can be computed using lattice QCD. The proposed approach requires well established 
methods, successfully applied in the calculations of electromagnetic and semileptonic form 
factors. The extra power divergences, related to the use of weak four-fermion operators, 
can be eliminated using only the symmetries of the lattice action without ambiguities or 
complicated non-perturbative subtractions. We demonstrate that this is true even when 
a lattice action with explicit chiral symmetry breaking is employed. Our study opens 
the possibility of reducing the present uncertainty in the theoretical predictions for these 
decays. 



1 Introduction 



Rare decays mediated by flavour-changing neutral-currents (FCNC) are among the deep- 
est probes to uncover the fundamental mechanism of quark flavour mixing. Within the 
Standard Model (SM), these rare decays are strongly suppressed both by the GIM mech- 
anism pP and by the hierarchy of the CKM matrix and are often dominated by 
short-distance dynamics. As a result, FCNC processes are very sensitive to possible new 
sources of flavour mixing, even if these occur well above the electroweak scale. The sensi- 
tivity to physics beyond the SM of these rare processes is closely related to the theoretical 
accuracy on which we are able to compute their amplitudes within the SM. 

Within the family of FCNC decays, long-distance effects are not always negligible 
and, in most cases, they represent the dominant source of theoretical uncertainty. Long- 
distance contributions are typically relevant in: i) amplitudes where the GIM mechanism is 
only logarithmic; ii) amplitudes where the power-like GIM suppression of the long-distance 
component is partially compensated by a large CKM coefficient. So far, the evaluation of 
these non-perturbative contributions has been performed by means of effective theories. 
These analytic tools require the introduction of additional parameters, the knowledge of 
which constitutes a source of sizable theoretical uncertainty. 

In this paper we show that for a class of very interesting processes, such as tt^uu 
and K — TTi~^i~, it is possible in principle to compute non-perturbatively the long- 
distance contribution to the physical amplitudes on the lattice. The physical information 
is encoded in the following T-products: 

rlJ{q')=^^v J d'x J d^'ye-^'^-y {n\T[Qix)J^{y)]\K) , (1) 

where Q denotes a generic four-quark operator of the effective weak Hamiltonian, J'^ 
is either the electromagnetic or the weak neutral current, and A/y is an appropriate 
volume factor. If the invariant mass of the lepton pair (g^) is smaller than any physical 
hadronic threshold, the calculation proceeds as in the case of semileptonic form factors 
(see e.g. refs. jSlIll), and one obtains directly the relevant amplitude. When instead the 
leptonic invariant mass exceeds the pion threshold, the final state interaction induces 
problems similar to those encountered with non-leptonic kaon decays [H] . However, the 
knowledge of the amplitude for < is sufficient to determine the leading unknown 
effective couplings of these amplitudes within the framework of chiral perturbation theory 
(CHPT) [ZI^[0!- Therefore, the combination of lattice calculations and CHPT should allow 
to reach an unprecedented level of precision for these rare decays. 

When using a lattice action with explicit chiral symmetry breaking, such as Wilson, 
Clover or twisted mass fermions, further problems arise because of additional ultraviolet 
(power) divergences which may appear in the operator matrix elements or in the relevant 
T-products.^ We show that for the electromagnetic current, gauge invariance prevents 

^ Alternative formulations which guarantee chiral symmetry in the physical matrix elements, such as 
overlap fermions, do not have this problem jlO| . However, these formulations are not mature yet to be 
used for unquenched calculations of these complicated matrix elements, for quark masses close to the 
physical values. 



1 



the appearance of these divergences even if the most popular lattice actions are used. 
Consequently, when J'^ is the electromagnetic current, the T-products in eq. ((T)) are fi- 
nite provided that a renormalized weak effective Hamiltonian is used. The situation is 
slightly more complicated when J'^ is the weak neutral current. In this case, simple power 
counting, related to the behavior of the T-product at short distances, shows that both 
quadratic and linear divergences may appear. We show that the quadratic divergence, 
which is not a peculiarity of the lattice regularization, is canceled by the GIM mecha- 
nism. Concerning the linear divergence, which is present only if there is an explicit chiral 
symmetry breaking term in the lattice action, we demonstrate that it can be avoided by 
using the maximally twisted mass fermion action jllj . 

There is a further subtlety concerning the ambiguity in the renormalization of the 
effective weak Hamiltonian out of the chiral limit jT2]. In ref. ^2] it has been shown 
that this ambiguity does not affect the physical K ^ amplitudes, but is present in 
"non-physical" matrix elements, such as (irlQlK). This problem is present also in our 
case and implies an ambiguity in the T-products of eq. (0). By means of appropriate 
Ward Identities, we show that the physical amplitude, the extraction of which requires a 
specific spectral analysis discussed in the following, is instead free of ambiguities. 

The paper is organized as follows: in sect.Elwe recall the basic ingredients of radiative 
decays in the framework of the effective Hamiltonian approach. In sect. 01 we describe the 
strategy for computing the relevant amplitudes from the Euclidean Green functions and 
discuss the structure of the divergences in both cases: when they cancel because of gauge 
invariance and when it is necessary to get rid of them using GIM mechanism and twisted 
mass. In sect. |3]we show how to extract the physical amplitude in spite of the ambiguity 
of the renormalized effective Hamiltonian. The results are summarized in the conclusions. 



2 Effective Hamiltonian for K iri^i (z/z/) decays 

The dimension-six effective Hamiltonian relevant to evaluate s di~^i^{h'u) amplitudes 
at next-to-leading order accuracy, renormalized at a scale Mw ^ /i > ttLc, can be written 
as 



He 



ntf-' + nlfr + ^ E + % E ^^.« + h.c., (2) 

q=u,d,s,c q=u,c q'=d,s 



where Vij denote the elements of the CKM matrix. 
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is the usual \AS\ = 1 weak Hamiltonian, for which the Wilson coefficients are known at 
the NLO CHI, and 
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Here 

Qj? = g7"(l-75)g'^>7M(l-75K 
gf = g7'^[2f(l-75)-4Qsin2^^]g 

X _£^^(i_^5_4sin2^H/K] (5) 

are the charged-current and neutral-current effective interactions obtained by the integra- 
tion of the heavy W and Z fields, 

(6) 
(7) 
(8) 

are the leading FCNC operators, and 

Ql = S„7''(l - 75)g/3 q(3lf,{l - l5)da , (9) 

Ql = Sa7''(l -75)?^ 9/37m(1 -75)c?/3 (10) 

the leading four-quark operators. The four-quark operators originated by penguin con- 
tractions are denoted by Qi,,,q, whereas Q7 and Qg correspond to magnetic and chromo- 
magnetic operators, respectively (see e.g. ref. [IH]). 

Thanks to both the GIM mechanism and the unitarity of the CKM matrix, the con- 
tributions to the FCNC amplitudes can be unambiguously decomposed into two parts, 
the first one proportional to the CKM combination V*gVud, the second one proportional 
to V^*Vtd- Since |Vj*Vtd| -C iKTsV^d], the contribution proportional to V^lVtd is negligible 
but for cases where it is enhanced by the large top-quark mass (i.e. for amplitudes which 
exhibit a power-like GIM mechanism). In these cases, the amplitudes are completely dom- 
inated by short distances (top-quark loops) and can be evaluated in perturbation theory 
to an excellent degree of approximation. In this paper instead we are interested only in 
the long-distance components of the amplitudes, therefore we can safely work in the limit 
Vtd = 0. 

We can seemingly neglect the matrix elements of Qi,„q and Qs in the evaluation of 
K — > tt£~^£~ (z/z/) amplitudes: these matrix elements vanish at the tree level and the 
corresponding Wilson coefficients are substantially smaller than those of Qi2- In this 
approximation, we only have to consider the contributions of the leading FCNC operators 
in eqs. (jHl)^® and the non-trivial contractions of 2 with the electromagnetic current 
and the currents defined by Q^J? and Q^'". The K ^ tt matrix elements of the FCNC 
operators in eqs. ©^(0) can be extracted from data on the leading K13 modes using isospin 
symmetry ^3], or even computed directly on the lattice, with high accuracy, as recently 
shown in Concerning the contractions of Q12, those with a charged current receive 

^ In principle, in the K Tr£~^£~ case one should also consider the tree-level matrix element of 
the magnetic operator Qy — ms'sa^^{l — j5)dF^,^, which cannot be directly extracted from Kg^ data. 
However, within the SM the sniallness of the corresponding Wilson coefficient makes this contribution 
negligible for practical purposes. This matrix element can be computed on the lattice with standard 
techniques, as shown in 
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very small non-perturbative contributions (estimated to be below 1% at the amplitude 
level in the — tx^vv case and even smaller in all the other channels), which can be 
reliably estimated within CHPT [HlllSI- Thus the main problem are the contractions of 
(5^2 with a neutral current, as outlined in eq. ((H). 

So far, this problem has been addressed with the following two-step procedure: i) 
integrating out the charm as dynamical degree of freedom; ii) constructing the chiral real- 
ization of the corresponding effective Hamiltonian with light quarks only. This procedure 
suffers from two sources of theoretical errors: slow convergence of perturbation theory be- 
cause of the low renormalization scale of the effective Hamiltonian (/x < mc); uncertainties 
associated to the new low-energy couplings appearing in the effective theory. Both these 
sources of uncertainties are naturally reduced in the lattice approach, where the effective 
Hamiltonian is renormalized above the charm scale and the T-products are evaluated in 
full QCD. 

We now discuss separately electromagnetic and neutrino amplitudes in more detail. 



2.1 K TT^+r 

The main non-perturbative correlators relevant for these decays are those with the elec- 
tromagnetic current. In particular, the relevant T-product in Minkowski space is [ZllHj 

= j d'xe-''''-{^={p)\T{J^^{x)m^)-Qm]]\KKk)) . (H) 
= (12) 

q=u,c q=d,s 

for 2 = 1,2 and j = +,0. Thanks to gauge invariance we can write 

(^') = ^ mk + pr - {ml - ml)q^] . (13) 

The normalization of (jl3p is such that the 0(1) scale-independent low-energy couplings 
a_|_^o defined in |8] can be expressed as 



1 

7^' 



2N- 
sm &w 



(14) 



where /+ is the K ^ tt vector form factor and {A^+,Ai'o} = {1,2~-^/^} To a good 
approximation, the decay rates of the CP-conserving transitions — »• 7i'~^i~^i~ and 
Ks 7T^i~^i~ are proportional to the square of these effective couplings |Hj: 

B{K+ TT+e+e-) ^ 6.6 x 10"^ , B{Ks ^ 7r°e+e") ^ 10.4 al x 10"^ . (15) 

At present, we are not able to predict a+,o with sufficient accuracy: we simply fit their 0(1) 
values from the measured rates of the corresponding decay modes (an updated numerical 
analysis can be found in [Ej). Being completely dominated by long distance contributions. 
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these two CP-conserving processes would provide an excellent testing ground for the lattice 
technique. 

On the other hand, the calculation of ao from first principles would have a very in- 
teresting phenomenological application in the Kl —>■ 7r^£~^£~ case, which proceeds via 
a CP-violating amplitude: the calculation of ao would allow to determine in a model- 
independent way the sign of the interference between the (long-distance) indirect-CP- 
violating component of the amplitude and the interesting (short-distance) direct-CP- 
violating term J7j. This result would allow to perform a very precise test of direct- 
CP-violation in the kaon sector. 

2.2 K TTVV 

The power-like GIM mechanism of the leading electroweak amplitude, implies a severe 
suppression of long-distance effects in these modes. In the CP-violating channel, — *■ 
n^uP, long-distance contributions are negligible well below the 1% level [16j. However, 
this is not the case for the charged channel, n^uu, where the suppression of 

long-distance effects is partially compensated by a large CKM coefficient. The T-product 
which determines the size of non-perturbative effects in this mode is [THl IT^ 

{T^r,{q') = -^ J d'xe-^'^-{n^{p)\T{J^,{x)[Qm-Qm]}\KHk)) , (16) 

where = g7''(2T(l— 75)— 4Q sin^ Ow)(l is the neutral current defined by Q^*-^ in (0). Sep- 
arating the electromagnetic component, we can write {%^)^ = (^^)i~4sin^ 9w (^^)em) 
where 

= -^Jd'xe-^'^-{7r^{p)\T{J^{x)[Qm-Qm]}\KHk)), (17) 
J'l = E^^^'^^i-^s)?- $^g7'^(l-75)g. (18) 

q=u,c q=d,s 

Contrary to (7^"^)|^^, the structure of ('^^)^ is not protected by gauge invariance and we 
can decompose it as 

{r^r^iq') = ^ [4ie){k+pr + o{qn] , (19) 

where the normalization is such that the zf{q'^) are expected to be 0{1) P^. The value 
of these form factors at = is sufficient to control the long-distance contributions to 
the —>■ TT^vv amplitude down to the 1% level of precision 0. 

Charm and, more generally, long-distance contributions to the — > tt^uu amplitude, 
are usually parametrized in terms of a scale-independent coefficient Pc ^Hl- According to 
the decomposition (fTI?|) . this can be be written as 

= ^ [C.ztiO) + C.ztiO)] + /+(0)a} . (20) 
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The coefficient Pc expresses the relative weight of the subleading terms relative to the 
top-quark amplitude, which is the leading contribution and is precisely determined in 
perturbation theory ^H]- As can be noted, the non-perturbative parameters z^{0) appear 
in multiplied by a very small coefficient: m|-/M^/|\4s|^ ^ 0.015. Thus even a 
determination of these matrix elements at the 30-50% level from lattice QCD would be 
sufficient to reduce the overall error on the n^uu rate around or below the 1-2% 

level. 



3 T-products at short-distances on the lattice 

In this section we discuss the properties of the Euclidean Green functions necessary to 
extract the physical amplitudes defined in eqs. (fTT|) and (fTTj) in a numerical simulation. 
Since the ultraviolet behavior is quite different in the two cases, we discuss them sepa- 
rately, starting from the T-product which involves the electromagnetic current. In both 
cases, we assume that the operators of the effective weak Hamiltonian have been renor- 
malized, namely that all their physical matrix elements are finite as the lattice spacing 
goes to zero (a 0). The renormalization of the effective Hamiltonian is discussed in the 
next section. 

The starting point to extract the physical matrix elements is the following Euclidean 
Green function 

{%)';,{q^U,tK) = I d'x{^^{U,p)J^{0)[Qt{x)-Qt{x)]^i{tK,k)) , 

U>0, tK<0, (21) 

where the source (sink) for creating (annihilating) the pseudoscalar mesons at fixed space 
momentum are defined as 

HU,qd = J d'ze-'^^-'^.iti.z), (22) 

and $i(ti, z) is a suitable local operator with the quantum numbers of the pion or kaon, 
respectively. Note that, in order to simplify the notation and the comparison between 
continuum and lattice formulae, we use the symbol of integral also to indicate sums over 
the lattice sites. 

If not for the presence of the weak four-fermion operator, the calculation would proceed 
as for the standard weak and electromagnetic form factors, by studying the behavior of 
the Green functions at large and This would give the form factors computed at 

momentum transfer q = k—p and with energy transfer go = Er — E-k- Since Qi is summed 
over the whole lattice volume and hence it carries zero momentum, this general strategy 
remains valid also for the Green function in eq. (j2^- As explained in the previous section, 
in order to extract the relevant low energy couplings, we are interested only to study the 
correlation function for < m\. In this range no rescattering of intermediate states is 
possible and thus we do not have problems in relating the Minkowskian T product to the 
Euclidean one. 
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Figure 1: One-loop topology which can originate power-like singularities to the Green 
function \21\) for x ^ 0. The dotted line denotes the generic insertion of Q^''^, with 
possible Fierz re- arrangements. 

The additional problem which arises in this case is the possibility that the Green 
function itself diverges because of the short distance behavior when x — 0. By dimen- 
sional arguments, this divergence can at most be quadratic. At fixed lattice spacing a, 
this would imply potential contributions to the Green function of (9(l/a^). Fortunately 
this never happens, since the strongest divergence associated to the diagram in figure ^ 
is independent of the quark masses and is canceled by the GIM mechanism. However, 
this cancellation does not guarantee the absence of linear divergences, which are naturally 
present when using lattice actions which break explicitly chiral invariance. 

3.1 The electromagnetic current 

Even if the chirality of the fermion action is explicitly broken, we are still able to define 
a conserved vector current on the lattice, which we can identify with the electromagnetic 
one. For example, with Wilson fermions we have 

J^ = ^[q{x + fi)U''\x){r + -f^')q{x)-qix)W\x){r--f^')q{x + fi)] , (23) 

where is the link variable. With a conserved current, gauge invariance is strong 
enough to protect the Green functions from the appearance of both quadratic and linear 
divergences. This remains true even when the Wick contractions correspond to a vacuum 
polarization diagram of the type in figure^ where only one of the two currents is the lattice 
conserved one, and the other is a local vector current originating from the weak four- 
fermion operator. We have verified this argument by an explicit perturbative calculation 
using Wilson, Clover and twisted mass fermions. Since the results of this calculation 
(more precisely of the subdiagram in figure ^ could be useful for other applications, we 
give them below for the Wilson and Clover cases. 
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p + q 



Figure 2: Subdiagram of figure^ associated to the x — > singularity. 



The amplitude we have considered is 

74 



^AV) = y^Tr [r«(g;p + g)A(g)7^A(g + p) 

= 7-^(W'-PMPv){X(p'a2,mV) + L} , (24) 
[An) 

where T']}\q] P+q) is the vertex derived from eq. ()23|) and A(g) the fermion propagator 
Both in Wilson and Clover cases we can identify a universal infrared term, given by 



T(p^a^, rn^a?) = dx x{l — x) log [m^a^ + p^a^a 
Jo 



a'x{l-x)\ , (25) 

while the finite constant L depends from the details of the regularization. In the Wilson 
case we find 

1 / 2 2\ , /i r X r A -^cosgpsin^g^ - |cosgpCos2g^ 
Lw = -g log (m a ) + (1 - 6,.) J_^^^ + 

(I cos Qp COS si^^ ^ ~ I COS 9'(T sin^ Q'p - | cos gp sin^ g^) 



^2 a 



^ A|(g;m) 
while in the Clover case 



r 



2 /•'T 1 <=;ri2 



sin^ gp - cos g^ sin^ ( y) 



^c.-^w + ^(l-MXA -""- ;7(;- (26) 

Note that in both cases the absence of power divergences holds independently from the 
GIM mechanism. Using the above results we could then match the lattice calculation 
with the continuum one even in an effective theory where the charm quark is integrated 
out. The comparison of the results obtained with or without dynamical charm quarks 
would provide a useful insight about the validity of the standard effective theory obtained 
by renormalizing Tie// below the charm mass. On the other hand, when the calculation 
is performed with a dynamical charm, the logarithmic divergence (and even the finite 
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coefficient) in eq. ()24|) is cancelled by the GIM mechanism. For this reason no matching 
lattice to continuum is needed in this case. 

Beside the possible singularities for s — > 0, further divergences may arise from contact 
terms of Qi with the external sources, namely for x x-,^ oi x ^ xk- However, it is easy 
to show that these contact terms do not contribute to the physical amplitudes. Let us 
consider the Minkowski T-product 

{%nxiQ'^i^^tK) = -t J d'x {0\t[^^{U,p)J^{0) [Q^{x) - Q^,{x)]^K{tK,k)] \0) , 

(27) 

corresponding to the Euclidean Green function of eq. ()21|) . The contact terms are propor- 
tional to the following pole terms: (p^ — m^)~^(A;^ — Tn'^)~^ or (p^ — m^)^^(A;^ — m^-)^^, 
while the on-shell amplitudes are obtained form the coefficient of {p^ —m'^)~^{k'^ — m\)~^ . 
As we shall discuss in more detail in the next section, these different pole structures in 
the Minkowski space correspond to a different — oo and tx ^ —oo behavior in the 
Euclidean case. As a result, we can eliminate the contact terms by an appropriate spectral 
analysis of the Green function computed in the numerical simulation. 

3.2 The axial current 

With the axial current appearing in the T-product (fTTj) . which is relevant for K — > -kvv 
decays, we cannot invoke gauge invariance: it remains true that the quadratic divergence is 
canceled by GIM, but we must face the problem of the linear one. With power divergences, 
any subtraction procedure, though non-perturbative, would produce an irreducible (and 
thus unacceptable) ambiguity in the final result. This implies that the linear divergence 
can only be an artifact of the regularization procedure. This divergence is indeed absent 
in regularizations which preserve chirality. 

With Wilson fermions the explicit breaking of chiral invariance leads to the appearence 
of such linear divergence. Since this problem is associated only to the contact term of the 
integrand PT|) for x ^ 0, we can in principle obtain a finite subtracted T-product, with the 
correct chiral behaviour of the Green function, by an integration which avoids the region 
close to X = 0.^ Otherwise, one could introduce an appropriate set of counterterms and 
fix their values by imposing an appropriate set of Ward identities, to recover the correct 
chiral behaviour. However, both these procedures are technically very complicated to be 
implemented. 

A much simpler and techincally feasible solution is obtained by means of maximally 
twisted mass terms Jl]. In this case, the additional symmetries of the action imply that 
the amplitude we are interested in is even in the Wilson parameter (r). This, in turn, 
implies the absence of the linear divergence which can only be odd in r, being associated 
to the breaking of chirality. We have verified this statement by an explicit perturbative 
calculation at the one loop level. As expected, the structure of the divergent terms is the 
same as in the continuum and the result is free from ambiguities. The discussion of the 

^We thank Massimo Testa for discussions on this point. 
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axial current can be repeated for a "non-conserved" vector current, such as the lattice 
local electromagnetic current, or the vector component of the weak left-handed current 
in eq. (fTTjl . 

At this point we wish to comment about the possibility to determine the physical 
K ^ TTTT amplitudes by using information about the following T-product 



where the axial current has the quantum numbers of the pion, but the kinematical 
configuration does not correspond to an on-shell pion field [221 • The Wick contractions for 
this T-product are similar to those considered for K —>■ nuu. In particular, the quadratic 
divergence generated when x — *■ is present also in this case. However, the situation is 
worse than the case discussed in this work, since there is no GIM mechanism to cancel 
the leading singularity. Of course we can define a renormalized T*-product, but this 
would entail a finite ambiguity. The practical problems which need to be faced in order 
to avoid this ambiguity make this calculation very difficult (if not practically impossible) 
with Wilson-type fermions. For this reason, we do not believe that a lattice study of 
this T-product can provide a useful tool to simplify the problem of determining K ^ titt 
amplitudes. 



In this section we address the problems arising by the renormalization of the lattice op- 
erators of the effective weak Hamiltonian. We first note that only the parity-even or 
parity-odd terms of the operators contribute to the vector or axial-vector cases, respec- 
tively. This observation is relevant since parity-even and parity-odd parts of the operators 
renormalize in a different way under regularizations which break chiral symmetry. On gen- 
eral grounds, whether chirality is broken or not, the mixing with operators of dimension 
five or six, in the presence of the GIM mechanism, does not introduce any ambiguity and 
the corresponding mixing coefficients can be computed in lattice perturbation theory. The 
problem arises from the mixing of the standard dimension-six operators with the scalar 
and pseudoscalar densities, which we now consider separately for the two cases. 
Schematically, we can write the renormalized operator as 

= Z^{^a)[Q^ + Cp{mc-mu){ms-md)s'y5d + Cs{mc-rnu)sd] , (29) 

where = (Qi ± (52)/2 + . . . represents the ensemble of all the dimension six and 
five operators with mixing coefficients computed in perturbation theory. By dimensional 
arguments, it follows that the coefficients Cs and Cp are power divergent in the limit 
a 0: 




(28) 



4 ^eff ambiguities 



(Jp — 

a 




(30) 
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Using suitable Ward identities (subtraction conditions), we can cancel the divergent parts 
of Cp and Cs', however, this leaves an ambiguity in their finite values out of the chiral 
limit p^l21j. For physical K — > vrvr amplitudes this ambiguity turns out to be irrelevant: 
the pseudoscalar density is proportional to the four divergence of the axial current and 
its matrix element vanishes for the on-shell K ^ transition jT2]- We stress that this 
conclusion does not hold for the K ^ n case: the off-shell matrix element {7T\sd\K) is 
different from zero thus, in general, the {7i\Q\K) matrix element does suffer from this 
ambiguity. 

Since we are interested in physical amplitudes, we must be able to demonstrate that 
also in the case of radiative decays the matrix elements of scalar and pseudoscalar densities 
do not contribute to the on-shell amplitudes. This can be done by means of suitable Ward 
identities and the spectral analysis of the relevant Euclidean Green functions. 

In the vector case we can use the following Ward identity 

d'^x [v^V;'^(x) + {nid - m,)srf(x)] J^(i/)<l>J^(xi^)) 

= -($i,(x.)J^(2/)$J,(x;,)) + ($^(a;^)J(;(l/)$t(xi,)), (31) 
where the term between square bracket is the rotation of the lattice action, 

= -^[s(x)f/^(x)(r-7^)f/(x + ^)-s(x + /x)f/;(x)(r + 7^)rf(x)], (32) 

and the last two terms in (j3T| correspond to the rotation of the pion sink and the kaon 
source, respectively. The term with the four divergence of V^'^(x), integrated over all 
space, vanishes. Thus on the left-hand side we are left with the term we are looking 
for, up to overall factors, namely the contribution of the scalar density to the Euclidean 
Green function (jH]), which enters when the bare weak operators are replaced with the 
renormalized ones. We need to show that this term does not contribute to the physical 
amplitude. 

In the Minkowski space, the physical amplitude is identified by the coefficient of the 
physical pole for —>■ and k"^ —>■ m|^. In the Euclidean space, this corresponds to a 
well-defined dependence on tx and (for ^ oo and tx — oo), namely 

^ ^ ^-EK\tK\ X e"-^-*- . (33) 



^p2 _ 7722)(A;2 _ mj^) 



The Ward identity (jHTj) tell us that the contribution of the scalar density give rise to a 
different pole structure: 



(p2 — m^)(A;2 — m^) 
1 

(p2 — m\-){k? — mj^) 



^-EK\tK\ X e"-^^*- (34) 
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Thus the scalar density contribution can simply be eliminated by a study of the time 
dependence of the appropriate Green function. Incidentally, this procedure eliminates 
also the divergent contact terms mentioned at the end of the section 13.11 

In the axial case, we have a similar situation, up to terms which vanish (linearly or 
quadratically in the lattice spacing) and inessential numerical factors. In particular, we 
can use the following Ward identity 

j d^x [-V^ZaA'^^{x) + (m, + ms)sj5dix) + 0(a)] rM^^j,{xK))] = 

= + J^(i/)St(x;,)) , (35) 

where again the term between square bracket is the rotation of the lattice action (for the 
explicit expressions of Za and the weak renormalized axial current see ^5) and Sj is a 
scalar particle source. This immediately shows that also the pseudoscalar density give 
rise to a time dependence different from the one in eq. and thus does not contribute 
to the on-shell amplitude. 



5 Conclusions 

The potential of rare K decays in performing precise tests of the SM and setting stringent 
bounds on physics beyond the SM depends, to a large extent, from our ability compute 
their amplitudes within the SM. In this paper we have shown that for a class of very inter- 
esting processes, such as -k^vv and K ni^i' , the theoretical error associated to 
non-perturbative effects could be reduced by means of lattice calculations. In particular, 
the numerical study of the Euclidean Green functions in eq. ()21|) . combined with CHPT, 
should allow to reach an unprecedented level of precision for these rare decays. 

The main problem which needs to be addressed before starting a lattice calculation 
of these Euclidean Green functions is the absence of power divergences in the extraction 
of the physical amplitudes. These may originate from contact terms between the weak 
four-fermion operators and the external fields (tt, K and the lepton current), or from the 
mixing of the four fermion operators with operators of lower dimensionality. In this paper 
we have shown that both these problems can be solved. 

As demonstrated in section HJ the spectral analysis necessary to extract the physical 
amplitudes eliminates both the power divergences due to the operator mixing and the 
contact terms with the external vr and K fields. The only remaining issue is then the 
ultraviolet behavior associated to the contact terms between the weak operators and the 
lepton current. This point is different for weak and electromagnetic currents. 

In the electromagnetic case, relevant for K —>■ TT£~^i~ decays, gauge invariance prevents 
the appearance of power divergences for all the popular Wilson-type actions. The can- 
cellation of power divergences is also independent of the GIM mechanism. We can thus 
match the lattice calculation with the continuum one also in an effective theory where the 
charm quark is integrated out. The perturbative expressions necessary for this matching 
at the one-loop level have been presented both for Wilson and Clover fermions. The 
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situation is slightly more complicated for the weak (axial or vector) current, relevant for 
TT^uu decays, where we cannot invoke anymore gauge invariance. One can cancel 
power divergences also in this case with Wilson-type fermions, but only using maximally 
twisted mass terms and taking advantage of the GIM mechanism. 

In summary, our analysis shows that the numerical study of the Green functions rele- 
vant for K — > 7ci~^i~ decays can be performed with any Wilson-type action, independently 
of the GIM mechanism. On the other hand, the study of ti^vv decays on the lat- 

tice requires a more sophisticated action: with Wilson-type fermions the only possibility 
is to use maximally twisted mass terms. We believe that these results opens a new field 
of interesting physical applications to the lattice community. 
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